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The  effects  of  the  diaphragm  shape  of  a  loudspeaker  on  the 
radiation  and  diffraction  characterisitcs  are  discussed  using  convex  and  con¬ 
cave  domes  in  an  infinite  baffle.  The  least  square  error  method,  one  of  the 
weighted  residual  methods,  is  employed  as  the  mathematical  tool  for  the 
solution  of  the  problems. 

The  results  show  that  these  characteristics  are  highly  dependent  on 
the  shape  of  the  diaphragm,  even  if  it  vibrate?  like  a  piston.  The  response  - 
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o£  the  concave  dome  has  a  wide  peak  due  to  the  cavity  resonance,  resulting 
in  higher  radiation  efficiency.  The  convex  dome  has  lower  on-axis  pressure 
response  in  the  same  region  due  to  the  dispersion  of  energy  to  the  off- 
axis  direction. 

The  diffraction  of  sound  from  a  concrete  ring  source  by  the  convex  and 
concave  domes  is  also  investigated  in  order  to  discuss  the  interaction 
between  the  loudspeaker  units  of  a  complete  system.  The  convex  dome  has 
much  larger  diffraction  effects  than  the  concave  dome,  especially  in  the 
high-frequency  region.  It  is  shown  that  even  the  concave  dome  produces  a 
discernible  amount  of  amplitude  distortion  of  sound  radiated  from  an 
adjacent  source. 

The  radiation  and  diffraction  phenomena  are  well  explained  by  the  use 
of  a  graphical  representation  of  the  energy  flow  and  pressure  distribution. 
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ABSTRACT 

The  effects  of  the  diaphragm  shape  of  a  loudspeaker  on 
the  radiation  and  diffraction  characteristics  are  discussed 
using  convex  and  concave  domes  in  an  infinite  baffle.  The 
least  square  error  method,  one  of  the  weighted  residual 
methods,  is  employed  as  the  mathematical  tool  for  the 
solution  of  the  problems. 

The  results  show  that  these  characteristics  are  highly 
dependent  on  the  shape  of  the  diaphragm,  even  if  it  vibrates 
like  a  piston.  The  response  of  the  concave  dome  has  a  wide 
peak  due  to  the  cavity  resonance,  resulting  in  higher 
radiation  efficiency.  The  con*'  x  dome  has  lower  on-axis 
pressure  response  in  the  same  region  due  to  the  dispersion 
of  energy  to  the  off-axis  direction. 

The  diffraction  of  sound  from  a  concentric  ring  source 
by  the  convex  and  concave  domes  is  also  investigated  in 
order  to  discuss  the  interaction  between  the  loudspeaker 
units  of  a  complete  system.  The  convex  dome  has  much  larger 
diffraction  effects  than  the  concave  dome,  especially  in  the 
high-frequency  region.  It  is  shown  that  even  the  concave 
done  produces  a  discernible  amount  of  amplitude  distortion 
of  sound  radiated  from  an  adjacent  source. 

The  radiation  and  diffraction  phenomena  are  well 

/ 

explained  by  the  use  of  a  graphical  representation  of  the 
energy  flow  and  pressure  distribution. 
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CHAPTER  I 


INTRODUCTION 


1.1  Background 

The  cone-type  loudspeaker  has  been  used  since  the  birth 

of  the  direct  radiator  loudspeaker.  The  reason  for  this  is 

the  mechanical  strength  against  the  force  given  through  the 

voice-coil  due  to  its  diaphragm  shape.  A  paper  cone, 

because  of  its  satisfactory  properties  such  as  relatively 

large  Young's  modulus-to-de ns i ty  ratio  (E/p^)  and  high  loss 

factor,  has  been  used  as  the  main  loudspeaker  diaphragm  of 

the  direct  radiator  loudspeaker.  But  as  far  as  paper  being 

used  as  a  diaphragm  material.  E/p  •  does  not  have  enough 

m  ; 

I 

margin  for  us  to  try  a  flat  or  a  very  shallow  diaphragm. 

The  situation  has  been  changed,  however,  since 

i 

composite  materials  with  large  E/p  i  have  come  into  use. 

m  | 

With  these  materials,  it  has  become  possible  to  try  any  kind 
of  diaphragm  shape,  even  the  flat  diaphragm.  At  this  point, 
we  naturally  come  to  the  question,  "What  is  the  best  shape 
for  a  radiator?"  The  diaphragm  shape  has  effects  on  both 
the  vibration  characteristics  and  the  radiation 
characteristics.  These  latter  effects,  however,  have  not 
had  much  attention  paid  to  them  up  until  now.  This  is  our 


motivation  to  investigate  the  relation  between  the  convex  or 
concave  diaphragm  in  an  infinite  baffle  and  its  radiation 
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character! a  tic s  . 

Other  interesting  phenomena  are  found  in  the 
interactions  between  loudspeaker  units  in  a  loudspeaker 
system.  The  sound  emitted  from  a  high-frequency  unit  is 
diffracted  by  a  (normally  concave)  low-frequency  unit.  This 
effect  causes  the  modification  of  the  on-  and  off-axis 
pressure  responses  of  the  high-frequency  unit.  This 
diffraction  is  also  responsible  for  an  amplitude-modulated 
distortion  of  the  hige-f requenc.y  sound  when  the  diaphragm  of 
the  low-frequency  unit  vibrates  with  a  large  excursion. 
These  static  and  dynamic  interactions  between  each  unit  of 
the  loudspeaker  system  are  important  topics  that  need 
investigation  if  improvements  on  loudspeaker  design  are  to 
he  made.  Another  example  of  this  kind  of  phenomenon  is  the 
diffraction  of  noise  by  a  convex  or  a  concave  object  on  a 
plane.  The  diffraction  of  noise  will  change  the  pressure 
distribution  around  the  object,  and  if  the  noise  source  is 
close  to  this  object,  it  may  also  change  the  total  radiated 
power  from  the  noise  source.  Thus,  the  diffraction  of  sound 
from  a  source  adjacent  to  a  convex  or  a  concave  object  in  an 
infinite  baffle  are  of  great  interest. 


1.2  Previous  Studies 

The  radiation  problems  from  convex  and  concave 

objects  are  quite  different  from  each  other  in  the  way  they 


are  solved . 


The  radiation  from  a 


convex  object  in  an 
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infinite  baffle  is  equivalent  to  the  radiation  from  a 
radiator  with  a  symmetric  radietor  in  the  other  semi- 
infinite  space*  On  the  other  hand,  the  radiation  problem 
from  a  concave  dome  in  an  infinite  baffle  must  be  solved  in 
a  semi-infinite  space.  The  studies  on  the  radiation  from  a 
radiator  of  finite  size  in  an  Infinite  space  are  reviewed 
first. 

The  radiation  from  a  complete  sphere  or  spheroid  is  the 
one  to  which  the  method  of  separation  of  variables  is 
applicable.  The  solution  to  this  kind  of  problem  is  found 
in  many  books  such  as  [l]-[5].  One  of  the  interesting 
researches  in  this  field  was  achieved  by  Ikegaya  [6].  As 
one  of  his  problems,  he  investigated  the  radiation  from  an 
oblate  spheroid,  assuming  that  it  can  be  well  approximated 
by  a  portion  of  a  sphere  in  an  infinite  baffle.  The  results 
show  that  the  on-axis  response  is  highly  dependent  on  the 
height-to-radius  ratio.  These  results  are  compared  with  our 
results  later. 

Three  distinct  methods  of  integral  forms  have  been 
reported  to  obtain  numerical  solutions  of  the  problems  of 
acoustic  radiation  from  an  arbitrary  body,  the  simple  source 
method,  the  surface  Helmholtz  integral  equation,  and  the 
internal  Helmholtz  integral  equation.  The  simple  source 
method  was  employed  by  Chen  and  Schweikert  [7],  and 
McCormick  and  Baron  [8].  Chen  and  Schweikert,  assuming  a 
distribution  of  simple  sources  over  the  radiator  surface, 
gave  the  formulation  of  the  radiation  problem  from  a  shell 


4 


in  an  infinite  medium.  McCormick  and  Baron  utilised  dynamic 
influence  coefficients  from  the  vibration  in  vacuum  of  tha 
shell  and  a  potential  source  theory  approach  for  the  fluid 
to  discuss  the  radiation  from  a  cylindrical  shell  of  finite 
length.  Chertock  [9]  investigated  the  sound  radiation  from 
a  body  of  arbitrary  shape  using  the  surface  .elmholtz 
integral  equation.  The  third  method  was  used  by  Copley 
[10],  who  adopted  the  internal  Helmholtz  integral  equation 
which  results  when  the  field  point  lies  within  the  surface 
of  the  radiator.  As  Copley  [11]  and  Schenck  [12]  pointed 
out,  however,  the  first  two  methods  fail  to  provide  unique 
solutions  when  the  frequency  is  approximately  equal  to  any 
of  the  characteristic  frequencies.  Schenck,  showing  that 
the  third  method  is  subject  to  similar  difficulties  and  has 
undesirable  computational  characteristic,  proposed  a 
modified  Helmholtz  integral  equation  formulation  in  the  same 
paper. 

A  quite  different  approach  was  applied  on  the  radiation 
problem  from  a  finite  cylinder  by  Williams,  Parke,  Moran, 
and  Sherman  [13].  They  expressed  the  pressure  field  in 
terms  of  spherical  wave  functions  with  unknown  coefficients. 
Then,  using  the  least  square  error  method,  the  unknown 
coefficients  were  determined  from  the  boundary  condition. 
The  same  method  was  used  to  solve  a  radiation  problem  of  a 
piston  set  in  a  sphere  by  Thompson  and  Skudrzyk  [14]. 
Fenlon  [IS]  discussed  the  radiation  field  at  the  surface  of 
a  finite  cylinder,  applying  the  method  of  weighted  residuals 
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on  Che  internal  Helaholtz  integral  equation.  Hunt.  Knittel, 
and  Barach  [16]  employed  the  finite  elenent  method  for  the 
vibrational  characteristics  of  the  elastic  body  and  the 
acoustic  pressure  field  enclosing  the  radiator.  They  used 
the  analytical  method  to  obtain  the  boundary  conditions  for 
this  mathematical  model. 

To  the  author's  knowledge,  very  few  studies  on  the 
radiation  frn  a  concave  object  in  an  infinite  baffle  have 
been  discussed.  Ohie,  Suzuki,  and  Shindo  [17]  discussea  the 
radiation  from  a  concave  sound  source  in  an  infinite  baffle, 
where  the  space  inside  the  truncated  cone  was  approximated 
by  the  stairwise  combination  of  thin  cylinders,  and 
cylindrical  wave  functions  were  used  in  each  cylinder.  The 
same  kind  of  problem  was  investigated  by  Shindo,  Kyono, 
Yashima,  Yamabuchi,  and  Kagaw..  [  1 8  j  ising  the  finite  element 
method.  The  same  method  was  used  Hy  Sakai,  Kyono,  Morita, 
Yamabuchi,  and  Kagawa  [  i  -j  J  for  the  investigation  of 
radiation  from  a  horn  loudspeaker. 

Problems  of  acoustic  diffraction  or  scattering  by  a 
sphere,  an  infinite  cylinder,  or  a  spheroid  are  found  in 
many  articles  such  as  [20]  —  [25]  ,  in  which  the  analytical 
method  is  applicable.  Another  technique,  the  integral 
equation  approach,  has  been  used  by  Burke,  Miller,  Poggie, 
Pjerrou,  Maxum,  and  Meechara  [26]  for  the  scattering  by  an 
elastic  body.  Hunt,  Knittel,  Nichols,  and  Barach  [27]  used 
the  finite  element  approach  to  acoustic  scattering  from  an 
elastic  circular  plate.  Most  of  these  problems,  however, 
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deal  with  the  diffraction  or  scattering  of  a  plane  wave, 
which  are  different  from  the  problem  of  interest. 

As  it  was  shown  above,  the  finite  element  method  has 
been  used  for  both  problems  of  radiation  from  convex  and 
concave  radiators.  The  method  similar  to  the  one  used  in 
[13],  however,  will  be  employed  for  the  present  redlatlon 
and  diffraction  problems  since  it  is  quite  close  to  the 
analytical  method  except  for  the  way  the  unknown 
coefficients  are  de  te  rail!  ad  . 


1.3  Statement  of  the  Problem 
The  main  object  lies  on  the  general  discussion  of  the 
relation  between  the  shape  of  the  diaphragm  and  its 
radiation  and  diffraction  characteristics.  For  the 
genera1 lty  of  the  discussion,  it  is  desirable  to  represent 
the  diaphragm  shape  as  simply  as  possible.  A  portion  of  a 
sphere  (which  is  called  "dome"  hereafter)  ic  suitable  for 
this  purpioe  because  it  is  simply  represented  by  one 
parameter,  the  height-to-radius  ratio.  To  simplify  the 
problem,  an  infinite  baffle  are  assumed  for  che  radiation 
and  diffraction  problems.  The  radiation  and  diffraction 
charscteristlcs  such  as  on-axis  pressure  response,  radiation 
impedanci ,  phase  response,  directivity  pattern,  and  near¬ 
field  sound  distribution  is  discussed.  It  is  reasonable  to 
assume  axisyrametry  of  the  sound  field  for  the  radiation 
problem  of  a  loudspeaker.  A  point  source  located  close  to  a 
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convex  or  a  concave  done  nay  be  a  more  suitable  node.!  for  a 
three-dimensional  diffraction  problem.  Initially!  however, 
only  an  axisymaetric  sound  field  is  treated.  Aa  a  sound 
source  in  the  diffraction  problem,  an  axisymaetric  circular 
line  source  (which  is  called  a  "ring  source"  hereafter)  is 
assumed.  This  makes  the  analysis  much  simpler  and  still 
gives  a  fundamental  understanding  of  the  problem. 

The  study  addresses  the  following:  (1)  for  the 
radiation  problem,  the  effect  of  the  height-to-radius  ratio 
of  the  convex  or  the  concave  dome  in  an  infinite  baffle  on 
its  radiation  characteristics;  and  (2)  for  the  diffraction 
problem,  the  effect  of  the  height-to-radius  ratio  of  the 
convex  or  the  concave  dome  in  the  infinite  baffle  on  the 
radiation  characteristics  of  the  axisymmetric  ring  source. 


CHAPTER  II 


RADIATION  FROM  CONVEX  AND  CONCAVE  DOMES 


2 . 1  Geoae t  r lea  1  Model 


2.1.1  Convex  Dome 

The  first  problem  is  to  investigate  the  radiation 
phenomenon  from  a  convex  dome  in  an  infinite  baffle  in  the 
semi-infinite  fluid  medium  with  density  p  and  velocity  c. 
The  existence  of  the  infinite  baffle  is  expressed 
mathematically  as  the  boundary  condition  that  the  normal 
velocity  on  the  baffle  surface  is  equal  to  zero.  The  same 
boundary  condition  is  achieved  by  assuming  another  symmetric 
dome  on  the  other  side  and  removing  the  infinite  baffle. 
For  the  convenience  of  mathematical  treatment,  this  radiator 
composed  of  these  two  domes  in  an  infinite  space  will  be 
used  as  the  radiator  model. 

The  crot  s  section  of  this  model  is  shown  in  Figure 
2.1.  The  center  of  the  radiator  is  the  origin  of  the 
rectangular  coordinate  system  (x,y,z).  The  rai'iator  is 
symmetric  about  xy-plane  and  axisyrametric  about  the  z-axis. 
Both  domes  are  described  as  a  portion  of  a  sphere  with 
radius  R.  The  sphere  of  the  right  side  dome  has  the  origin 
0).  The  polar  coordinate  systems  (r'.d',^')  with  the  origin  0 
and  (r,6,4>)  or  with  the  origin  0^  are  used. 


The 
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area  of  the  right  side  dome  (r«R,  O<0<0q,  0<4><2tt)  is  denoted  by 

S^\  and  the  imaginary  surface  (r»R,  0  <0<t,  0<4»<2 tt)  is 

o  - 

denoted  by  The  radius  and  the  height  of  the  dome  are 

A  and  H,  respectively.  The  domes  are  assumed  to  vibrate 
sinusoidally  with  angular  frequency  ta  in  the  z  direction  and 
with  the  amplitude  UQ  and  180  degrees  out-of-phase  with  each 
other.  The  assumption  cf  the  piston-like  motion  may  not  be 
suitable  for  an  actual  loudspeaker  in  the  high-frequency 
region  because  of  the  higher  order  modes  developing  on  the 
diaphragm,  but  it  will  be  beneficial  for  basic  discussions 
of  sound  radiation. 

2.1.2  Concave  Dome 

Figure  2.2  shows  the  geometry  of  the  concave  dome 

radiator.  The  xy-plane  coincides  with  the  baffle  surface, 

and  the  z-axis  is  the  axis  of  revolution.  The  dome  is 

described  as  a  portion  of  a  sphere  of  radius  R  and  origin 

Oj.  The  dome  has  the  radius  A  and  height  H.  The  polar 

coordinate  system  (r.O.if)  will  be  used  for  the  semi-infinite 

space  including  the  cavity  in  the  baffle.  The  surface  of 

the  dome  (r*R,  9  <9<  tt,  0<c+)<2n )  is  denoted  by  The  points 

a~  — 

on  the  surface  of  the  opening  S  ^  ^  will  be  described  by 


(r1,01,<{)1)  . 
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2.2  Ma  t hema t ica 1  Discussions 


2.2.1  Basic  Equations  f or  Radiation  Problem 

The  sound  in  a  fluid  medium  is  a  dynamic  disturbance 
of  the  fluid  which  can  be  described  by  a  time  dependent 
scalar  quantity  like  pressure  or  velocity  potential.  In 
linear  acoustics,  the  wave  equation  of  the  harmonic  wave 
propagation  in  the  fluid  is  expressed  by  the  Helmholtz 
equa  t ion : 

(V2  +  k2)  ttf)  -  0  ,  (2.1) 

where  i|/(£)  Is  the  velocity  potential,  k  is  the  wavenumber, 

V2  is  the  Laplacian  operator,  £  is  the  field  point  in  the 
acoustic  medium,  and  the  time-dependent  term  e^1*11'  is 
suppressed.  The  pressure  p(£)  and  the  particle  velocity 
u(£)  are  related  to  the  velocity  potential  by  the  following 
equations : 

p(£)  -  jwp  i K£)  (2.2) 

and 

u(£)  -  -  V  i|>(£)  ,  (2.3) 

/here  V  is  the  gradient  operator. 

The  radiation  problem  is  to  find  a  solution  to 
Eq .  (2.1)  to  satisfy  the  appropriate  boundary  condition. 
The  boundary  condition  on  the  radiator  surface  Is  expressed 


by 
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li'te)  «  f(C)  ,  (2.4) 

where  5  is  a  point  on  the  radiator  surface.  The  Dirichlet 
and  Neumann  boundary  conditions  are  given  by 

L  *  jup  ;  f(0  -  ?(0  (2.5) 

and 

L  -  -  ~  l  f(t>  "  U(0  ,  (2.6) 

3n 

respectively,  where  P (?)  and  U(£)  are  the  pressure  and  the 
normal  velocity  distributions  prescribed  on  the  surface  of 
the  radiator,  and  IT  is  the  unit  outward  normal  to  the 
surface.  It  should  be  noted  that  Eq  .  (2.4)  can  also  express 

the  mixed  boundary  condition  if  one  assumes  Dirichlet  and 
Neumann  boundary  conditions  on  the  different  areas  of  the 
radiator  surface. 

The  solution  of  Eq .  (2.1)  also  has  to  satisfy  the 
radiation  condition  [28]: 


lira 

A  _ 

^3SOO 


9(0  +  jk  p(C) 


2 

r=>r 


dS 


0  , 


(2.7) 


r 

where  r  is  the  distance  from  the  origin  of  coordinates  to 
the  field  point,  and  S^,  is  the  surface  of  a  sphere  of  radius 
f  centered  at  the  origin,  surrounding  the  field  point  £  and 


the  radiator 
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2.2.2  Least  Square  Error  Method 

As  it  was  mentioned  before,  one  cannot  apply  the  method 
of  separation  of  variables,  because  there  exists  no 
appropriate  coordinate  system  that  matches  the  boundary  of 
the  radiator  except  in  the  case  of  H/A  ”  1.0,  which 
corresponds  to  the  radiation  from  a  complete  sphere.  Here, 
the  method  called  "least  square  error  method"  will  be  used 
extensively  in  the  present  radiation  and  diffraction 
problems  [ 29 ] . 

The  velocity  potential  will  be  expanded  in  the 
following  form: 


HO  "  Z  an  *  (O  ,  (2.8) 

n=0 

where  an  is  the  unknown  coefficient,  and  the  set  of  infinite 

crial  functions  ^  ( O  is  considered  to  be  capable  of 

n 

describing  the  unknown  function  which  satisfies  the 

wave  equation  Eq.  (2.1)  and  the  radiation  condition 
Eq  .  (2.7).  When  the  maximum  of  the  order  n  is  truncated  by 
N,  the  "best"  approximation  of  the  boundary  condition 
Eq .  (2.4)  on  the  boundary  S  is  defined  as  a  solution  with 
the  set  of  unknown  coefficients  an,  n  =  0,  1,  .  .  N  that 

makes  the  following  functional  stationary: 


r 

N 

J  = 

J 

q(C) 

S 

Z  a  w  (O  -  f(£) 

n=0  n  n 

(2.9) 


where 
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V°  "  L  ’  (2.10) 

a’  t:l  q(C)  is  a  non-negative  weighting  function  necessary  at 
least  for  the  consistency  of  dimensions  of  the  integrands  of 
Eq .  (2.9)  in  the  case  of  mixed  boundary  condition  problems. 
The  an  is  given  by  the  condition  that  J  is  stationary,  i.e., 
insensitive  to  arbitrary  variations  in  the  parameter  an, 


(2.11) 


Substituting  Eq.  (9)  into  Eq.  (11),  (N+l)  simultaneous 

equations  are  obtained  [see  Appendix  A]: 


S 


a  t  (0^(0 
n  n  m 


dS 


q(Cn'*(C)f(0dS 


m  =  0,  1,  .  .  . ,  N  , 


(2. 12) 


where  *  denotes  the  complex  conjugate. 

Thus,  if  the  proper  trial  functions  and  expressions 
for  the  boundary  condition  are  found,  the  solution  of  best 
approximation  can  be  obtained  by  solving  Eq  .  (2.12). 


2.2.3  Radiation  from  a  Convex  Dome 

The  general  discussion  of  the  mathematical  tool  used 
for  the  radiation  problem  has  been  completed  in  the  previous 
two  sections.  In  order  to  solve  the  problem  efficiently,  it 
is  important  to  use  physical  and  mathematical  intuition  for 
the  selection  of  appropriate  trial  functions  and  the 
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coordinate  system.  This  means  that  the  proper  choice  of 
trial  functions  and  coordinate  system  enables  one  to 
approximate  the  boundary  condition  accurately  enough  with 
the  smallest  order  of  truncation  N. 

The  direct  and  simple  application  of  the  previous 
mathematical  discussion  for  the  radiation  problem  by  the 
convex  dome  shown  in  Figure  2.1  is  to  expand  the  velocity 
potential  in  terms  of  spherical  harmonic  functions  with  the 
origin  at  the  center  of  the  radiator.  However,  two  basic 
problems  are  inherent  with  this  procedure.  [These  are 
discussed  in  detail  in  Appendix  B.) 

Now,  on  the  basis  of  this  knowledge,  it  seems 
reasonable  to  choose  the  center  of  one  of  the  spheres  as  the 
origin  of  the  coordinate  system.  The  velocity  potential  is 
expanded  in  terms  of  the  infinite  series  such  as  [30]: 


’Kx) 


CO 


(2.13) 


where  hn(kr)  is  the  spherical  Hankel  function  of  the  second 
kind  and  order  n  with  the  suppression  of  the  superscript 
(2),  and  P^CcosO)  is  the  Legendre  function  of  the  first 
kind  and  order  n.  [The  validity  of  the  expression  in  the 
form  of  Eq .  (2.13)  for  Dirichlet  and  Neumann  boundary 
conditions  is  shown  in  Appendix  C.] 

Then,  the  boundary  condition  on  the  radiator  surface 
S<1>  is  given  by 


(1) 


r~R 


9 

9r 


(2.14) 
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k 


r. 


and 

f(1)(c)  -  Uocose  .  I  (2.15) 

Then , 

’J,<1)(0)  -  -kh’(kR)  P  (cos6)  .  (2.16) 

n  n  n 

Here,  the  superscript  (1)  or  (2)  indicates  the  functions 

with  regard  to  the  surface  or  respectively. 

Also,  (6)  and  f^(9)  are  used  instead  of  4*^(0  and 

n  n 

f^(S)  since  they  are  functions  of  9  only. 

On  the  imaginary  surface  one  can  use  the 

condition  of  symmetry  of  the  sound  field.  The  velocity 
potential  at  point  P  is  equal  to  the  one  at  point  Pj,  giving 
the  relation: 

N  N 

l  a  h  (kr)  P  (cosQ)  =  l  a  h  (kr. )  P  (cos9  )  ,  (2.17) 

Oil  11  11  11  11  X  11  -1 

n=0 


where  R,  0,  r^,  and  0^  are  related  to  each  other  by  the 

following  equations: 

r^  cos0^  =  R  (2cos0q  -  cos0)  (2.18) 

and 


rx  sin01  =  R  sin0  ,  (2.19) 

and  R  is  given  by 

R  =  (H2  +  A2)/2H  .  (2.20) 


~K 
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Rewriting  Eq.  (2.17),  one  obtains 


N 

l  a  [h  (kR)  Pn(cos6)  -  h^kr^  Pn(cos61>]  -  0  .  (2.21) 

n-0  n  n  n 


The  condition  of  symmetry  of  the  sound  field  does  not 
give  an  explicit  form  to  the  operator  L.  But,  comparing 
Eq .  (2.21)  with  Eq.  (2.4),  one  comes  to  the  interpretation 
that  the  operator  L,  when  applied  on  the  velocity  potential 
at  a  point,  gives  the  difference  of  velocity  potentials  at 
that  point  and  its  symmetric  point. 

Thus,  and  f^(6)  can  be  defined  such  as: 

Hf^(0)  -  hn(kr)  Pn(cos6)  -  hn(kr1)  P^cosO^)  (2.22) 

and 

f(2)(0)  -  0  .  (2.23) 

Applying  Eqs.  (2.16),  (2.17),  and  Eqs.  (2.22),  (2.23)  into 

Eq .  (2.12),  the  set  of  (N+l)  simultaneous  equations  are 
obtained : 


N 

J  a  l 

n-0  n  > 


r  m  y(1)*(0>  ,^1)(0)  ds  + 
.(l)  m  n 


.(2) 


H'(2)*(0)  7(2)(0)  dS  1 
m  n 


'i'(1)*(0)  f(1)(6)  dS  ,  m-  0,  1,  ...»  N,  (2.24) 
m 


.(1) 


whe  re 


q(C)  = 


1.0 


(2.25) 
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Is  assumed  In  Eq .  (2.12). 


2.2.4  Radiation  from  a  Concave  PggLg. 

In  the  case  of  radiation  from  a  concave  dose,  it  is 
impossible  to  get  rid  of  the  infinite  baffle  without 
changing  the  boundary  condition.  Thus,  the  sound  field  must 
be  dealt  with  in  a  semi-in f ini t e  space.  Instead,  one  can 
use  the  knowledge  that  the  sound  field  is  uniquely 
determined  by  the  use  of  the  Rayleigh  integral  if  the  normal 
velocity  distribution  on  the  baffle  surface  is  given. 

The  velocity  potential  inside  the  sphere  with  the 
origin  Oj  is  expanded  in  terms  of  infinite  series  such  as 
[30]  : 


N 

YU)  ”  l  \  Jn(kr)  Pn(cos0)  » 

n-0 


(2.26) 


where  jn(kr)  is  the  spherical  Bessel  function  of  the  first 
kind  and  order  n- 

The  boundary  condition  on  the  radiator  surface 
is  given  by  the  same  forms  as  Eqs.  (2.14)  and  (2.15): 

T(l)  _  _  J_  (2.27) 

L  "  9r  r-R 

and 


f 


(1) 


(9) 


!J  cosQ  • 
o 


(2.28) 


Thus  , 
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«?(1^(0)  -  -  k  j'(kR)  P  (cosQ)  . 

n  n  n 


(2.29) 


Once  the  expression  for  the  velocity  potential  (with 

unknown  coefficients)  inside  the  sphere  is  given,  one  can 

(  2  ) 

obtain  the  normal  velocity  distribution  on  the  opening  S 
of  the  cavity: 


uB(rl‘0l>  “ 


3 


-£*(  r’9> 


r-r. 


e-9 : 


34,36]  . 

Dr  3z  39  3zj  1 


(2.30) 


9-6, 


The  derivatives  are  given  from  the  geometry  by 


3r 

3z 


COS0 


(2.31) 


and 


36 

3z 


sin9 

r 


(2.32) 


Substituting  Eqs.  (2.31)  and  (2.32) 
obtains 


into  Eq.  (2.30),  one 


wv 


l  a  [kcos91  jH krx)  Pn(cos01) 
n-0  n 


sin2©,  j  (kr.)  P*(cos9  )/r  ]  . 
I  n  1  n  l  J- 


(2.33) 


Now,  the  velocity  potential  ^(r.0)  due  to  ^(r^Oj)  outside 
be  obtained  using  the  Rayleigh  integral: 


.  s'x* 


the  cavity  can 
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1  f  f  e“jkd 

*u(r'e)  ■  2?  ,JU-<C1-61>  -3—  dS  •  (2.34) 

S(2)  n 

where 

2  2 

d  -  [(rsin0  -  r^sinO^cosifc^)  +  (r^sinO^sinii^) 

2  1/2 

+  (rcos0  -  r^cos0^)  ]  (2.35) 

is  the  distance  between  the  point  of  interest  (r,0,O)  and 
the  point  of  surface  element  dS  (r^.fl^,^)  .  The  two 
velocity  potentials  given  by  Eqs.  (2.26)  and  (2.34)  must  be 
equal  to  each  other  (potential  matching)  everywhere  inside 
the  sphere  and  outside  the  cavity. 

It  seems  reasonable  to  impose  the  condition  of  the 
potential  matching  on  the  imaginary  surface  of  the  sphere. 
But  the  expansion  of  the  form  of  Eq .  (2.26)  causes  a 
nonuniqueness  problem  when  the  boundary  condition  is  given 
on  the  surface  of  the  sphere  at  or  near  its  characteristic 
frequencies  [see  Appendix  C].  For  this  reason,  the  opening 
of  the  cavity  Sv<’/  is  chosen  as  the  area  on  which  the 
condition  of  potential  matching  is  imposed.  Then,  one 
obtains 

N  N 

l  an[jn(kr)  Pn(cos0)I  -  *  J  a  t  {kcose.V  (kr, ) 

n~0  n-0  n  J  (2)  1  n  1 

s 

2  -jkd 

X  IMcose^  +  sin  0X  3^)  P^ccs^)/^)  dS)  . 


(2.36) 
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rr 


with  the  relationship  between  r  and  6  such  as: 


r  ■  -  z  /cos0  . 

o 


(2.37) 


ih 

-k 

s 

ft 

*v! 


where  z  is  the  distance  between  0  and  0,.  Rewrite  this 
o  i 

equation,  "  - — 

N  f 

J  a  [j  (kr)  P  (cosO)  +  (kcosQ  V(kr  )  P  (cosQ.) 
n-0  n  n  n  Js(2)  1  n  x  n  1 

-jkd 

+  sinSj^  ja(kr1)  P^cosQ^/r^  (—p)  dS]  -  0  .  (2.38) 


Thus,  ¥^2y(9)  and  f^(9)  can  be  defined  such  as 


^(2)  (9)  -  [ j  (kr)  P  (cos9)  + 

n  n  n  j 


.(2) 


(kcosQ. j ’ (kr. )  P  (cos9,) 
l  n  x  n  x 


-jkd 


+  sin201jn(kr1)  P^(cos01)/r1)  (^f)  dS]  (2.39) 


and 


f(2)(6)  -  0 


(2. AO) 


Now,  by  the  same  method  used  in  Section  2.2.3,  one  obtains 
the  (N+l)  simultaneous  equations  [again  Eq  .  (2.25)  was 
assumed ] : 


J  a  [f  4,(1)*(9)  H'(1)(9)  dS  +  f  (9)H,(2) 

*•  n  L  I  m  n  '  '  IH  n 


n-0  n  Js(l)  m 


1  (2) 

SK  J 


(9)  dS] 


V1* 


<1>,V(°)f  (1>(°)  dS  *  «  -  0.  1,  .  .  .,  N  . 


(  2  .  A  1  ) 


23 


2.2.5  Geometrical  App rox Ima t Ion  Method 

It  is  well  known  that  sound  will  be  radiated  normal  to 
the  surface  of  the  radiator  when  the  wavelength  is  much 
smaller  than  the  radiator  size.  For  higher  frequencies,  the 
geometrical  approximation  begins  to  have  practical  meaning. 
The  geometrical  approximation  here  assumes  that  an  infinite 
number  of  point  sources  are  distributed  in  such  a  way  that 
they  form  the  shape  of  the  radiator  with  equal  surface 
density  in  an  infinite  space.  The  on-axis  sound  pressure  of 
the  convex  dome,  which  is  normalized  by  the  sound  pressure 
due  to  a  single  point  source  with  the  same  total  volume 
velocity  at  the  origin  0,  is  then  given  by 

•6  ikR(cos6-cos0  )  „  „ 

p  *  [  °  2iTRsin0cos0e  0  Rd9]/irR  sin  0 

Jq  ° 

ikR(l-cosQ  )  ?  ? 

-  2[e  °  (ikR  -  1)  +  (1  -  cos0Q  •  ikR)/sin  0Q(ikR)  . 

(2.42) 

The  expression  for  the  concave  dome  is  the  same  except  that 
the  phase  response  has  the  opposite  sign.  The  results 
obtained  with  this  method  will  be  compared  with  exact 


solutions . 
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2 . 3  Results  and  D1 scuss Ions 
The  calculations  were  made  for  three  different 
height/radius  ratios  such  that  H/A  *  0.5,  0.75,  and  1.0. 

The  frequency  range  was  from  kA  *  0.1  to  10.0.  The  accuracy 
of  the  results  was  confirmed  by  verifying  the  degree  of 
matching  of  velocity  and  velocity  potential.  Except  in  the 
vicinity  of  the  rim  of  the  dome,  the  differences  were  less 
than  a  few  percent.  The  weighting  function  q(C)  in 

Eq.  (2.12)  was  kept  equal  to  unity  all  through  the 
calculations,  because  varying  it  did  not  give  any 
significant  improvement  in  accuracy.  The  maximum  order  of  N 
was  40  for  kA  *  10.0. 

2.3.1  On-axis  Pressure  Response 

The  far-field  on-axis  pressure  response  of  a  convex 
dome  normalized  to  the  pressure  response  of  a  flat  piston 
with  same  volume  velocity  are  shown  in  Figure  2.3,  for  H/A 

■  0.5,  0.75,  and  1.0.  When  H/A  is  equal  to  0.5,  the 
response  shows  a  meaningful  decrease  from  about  kA  **  0.5, 
and  reaches  to  the  -3.5dB  line  around  kA  *  2.0.  For  H/A  » 
0.75  and  1.0,  minimum  levels  in  the  middle  frequency  region 
(say,  kA  ■  0.5  -  5.0)  are  -6.5dB  and  -9dB,  respectively. 
This  effect  of  the  diaphragm  shape  on  the  pressure  response 
is  essential  when  a  dome-type  loudspeaker  is  designed.  For 
a  dome-type  loudspeaker  with  a  1cm  radius,  the  range  from  kA 

■  0.5  to  5.0  corresponds  to  the  frequency  range  from  2.7kHz 
to  27kHz.  This  range  includes  the  reproduction  range  of  a 


Figure  -.3  Far-field  on-axis  pressure  responses  of  a  convex 
dome  radiator. 
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dome-type  loudspeaker,  which  may  have  the  H/A  ratio  of  about 
0.75.  This  means  that  the  achieved  pressure  response  is 
much  lower  than  the  one  estimated  from  the  assumption  of  a 
flat  piston.  This  kind  of  discrepancy  is  usually 
encountered  when  a  dome-tyoe  loudspeaker  is  designed. 

It  is  meaningful  to  compare  the  pressure  response  with 
results  obtained  by  the  geometrical  approximation  method, 
which  are  shown  in  Figure  2.4.  All  three  responses  begin  to 
roll  off  from  about  kA  ■  3.0  at  the  rate  of  -6dB/oct.,  but 
with  different  undulations.  As  the  H/A  ratio  becomes 
smaller,  the  fluctuation  becomes  larger.  In  Figures 
2. 5-2. 7,  two  responses  obtained  by  the  least  square  error 
method  (LSM)  and  the  geometrical  approximation  method  (GAM) 
are  compared  for  each  H/A  ratio.  Figures  show  that  the 
geometrical  approximation  gives  a  fairly  good  estimation  of 
the  true  response  in  the  high-frequency  region.  Especially 
for  H/A  *  1.0,  the  difference  is  less  than  ldB  above  kA  » 
4.0. 

The  results  of  the  concave  dome  are  shown  in  Figures 
2.8-2.10.  Contrary  to  the  results  of  the  convex  dome,  the 
pressure  response  of  the  concave  dome  has  a  wide  peak 
around  kA  *  1.0  -  1.5  for  H/A  ■  1.0  -  0.5.  This  region  is 
included  in  the  reproduction  range  of  an  actual  cone-type 
loudspeaker.  For  H/A  ■*  0.5,  the  peak  level  is  about  4dB. 
This  kind  of  resonance  effect  is  not  clearly  observed  in  the 
pressure  responses  of  an  actual  loudspeaker.  Two  reasons 
for  this  can  be  considered.  The  first  one  is  that  the 


approximation 
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actual  cone  diaphragm  may  break  up  around  this  frequency 
region.  The  second  one  is  that  the  inductance  of  the  voice- 
coil  reduces  the  driving  force  in  the  same  region.  For  the 
concave  dome,  the  geometrical  approximation  gives  a  rough 
estimation  of  the  true  response. 

The  above  discussion  shows  that  the  total  range  may  be 
divided  into  three  regions.  The  region  in  which  the 
response  is  almost  equal  to  that  of  a  flat  piston  is  the 
low-frequency  region  (approximately  kA  0.5).  In  the  high- 
frequency  region  (approximately  kA  >  5.0),  the  response  can 
be  estimated  by  the  geometrical  approximation  method.  In 
the  range  between  these  two  regions,  the  response  must  be 
obtained  by  an  exact  method. 

The  fact  that  responses  of  the  convex  and  concave 
domes  are  intrinsically  different  may  mean  that  the  sound 
quality  of  them  may  also  have  basic  differences. 

2.3.2  Radiation  Impedance 

The  normalized  radiation  impedance  is  defined  by  the 
following  equation: 

z  -  4  [  p(C)u*(C)dS}/{i  pcttaV)  .  (2.43) 

r  l  »s(1)  2  o 

The  numerator  of  this  equation  is  the  total  power  (complex), 

and  the  denominator  is  the  effective  radiated  power  from  a 

radiator  with  velocity  UQ  and  constant  radiation  impedance 
2 

pcxA  .  The  real  part  (Re)  of  zf  is  the  normalized 
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radiation  resistance,  and  the  imaginary  part  (Im)  is  the 
normalised  radiation  reactance. 

The  radiation  impedance  characteristics  of  the  convex 
dome  are  shown  in  Figures  2.11  and  2.12.  The  radiation 
resistance  of  the  convex  dome  is  smaller  than  that  of  the 
flat  piston  with  the  same  radius,  and  decreases  as  the 
height  of  the  dome  Increases.  The  radiation  resistance  of  a 
flat  piston  fluctuates  above  kA  ■  2.0,  while  its  pressure 
response  stays  constant.  On  the  other  hand,  for  the  same 
frequency  region,  the  radiation  resistance  of  the  convex 
dome  is  constant  even  if  its  on-axis  pressure  response 
varies  . 

In  order  to  discuss  the  frequency  response  of  the 
radiated  power,  the  radiation  impedances  were  plotted  on  a 
logarithmic  scale  in  Figures  2.13  and  2.14.  When  the 
diaphragm  is  driven  with  constant  acceleration  with  respect 
to  frequency,  the  velocity  decreases  by  6dB/oct.  If  the 
response  of  lQlog(*r)  increases  by  6dB/oct.,  the  radiated 
power  remains  constant.  As  Figure  2.13  shows,  the  flat 
piston  has  the  closest  radiation  resistance  to  the  line  with 
a  6dB/oct.  slope.  As  the  H/A  ratio  increases,  the  convex 
dome  reduces  its  range  of  the  constant  power  response.  For 
H/A  ■  0.5,  the  radiated  power  at  kA  ■  2.0  is  about  5dB  lower 
than  the  level  in  the  low-frequency  range.  The  radiation 
reactance  of  the  convex  dome  has  a  slope  of  +3dB/oct.  and 
-3dB/oct.  in  the  low-  and  high-frequency  region, 
respectively,  and  always  positive  (mass-like). 


onvex  done  and  a  flat  piston  (normalized  by 


-4tC3*UL • 


Figure  2.12  Radiat 
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The  radiation  impedances  of  the  concave  done  are  shown 
in  Figures  2.15-2.17.  Radiation  resistance  of  the  concave 
dome  is  much  larger  than  that  of  the  flat  piston  and  the 
convex  dome  in  the  region  kA  ^  3.0.  This  can  be  explained 

as  a  resonance  effect  of  the  cavity.  The  frequency  of  the 
maximum  resistance  is  slightly  higher  than  that  of  the 
maximum  pressure  response.  The  radiation  reactance  is  zero 
at  these  frequencies,  and  mostly  negative  (spring-like) 
above  those  frequencies. 

The  plot  of  the  radiation  resistance  on  a  logarithmic 

scale  is  shown  in  Figure  2.18.  The  6db/oct.  line  stays 

between  the  curves  of  the  flat  piston  and  concave  dome  with 

H/A  ■  0.3.  Thus,  if  the  acceleration  of  a  diaphragm  is 

constant  with  respect  to  frequency  and  if  the  diaphragm 

stays  like  a  piston,  the  widest  range  of  constant  radiated 

power  is  achieved  by  a  concave  dome  with  H/A  ratio  somewhere 

between  0.0  and  0.5.  The  radiation  mass  loaded  on  the 

concave  radiator  is  almost  equal  to  the  sum  of  the  mass  of 

the  fluid  inside  the  cavity  and  the  radiation  mass  of  a  flat 

3 

piston  with  the  same  radius  (*  0.85  pirA  ). 

2.3.3  Phase  Response 

The  far-field  on-axis  phase  responses  of  the  convex 
and  concave  domes,  which  are  normalized  to  the  phase 
response  of  an  on-axis  point  source  on  the  baffle,  are  shown 
in  Figures  2.19-2.20,  respectively.  The  phase  responses  of 
the  convex  and  concave  domes  increase  or  decrease 
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Figure  2.15  Radiation  impedance  characteristics  of  a 
concave  dome  with  H/A  =  0.5  (normalized  by 

pcnA^  }  ♦ 


Figure  2.17  Radiation  impedance  characteristics  of  a 
concave  dome  with  H/ A  =  1.0  (normalized  by 


4? 


respectively  as  the  frequency  increases.  It  is  interesting 
to  separate  the  total  phase  response  into  two  parts,  i.e., 
the  minimum  phase  response  obtained  from  the  amplitude 
response  through  a  Hilbert  transform  and  the  remaining  phase 
response.  The  transfer  function  of  the  loudspeaker  can 
therefore  be  represented  by  the  following  form: 


H(w) 


id)  (to)  id)  , ,  (til) 
.  ,  s  vminv  ’allpass 
A (w)e  e 


(  2  .  A  A  ) 


The  non-minimum  phase  response  ^allpass^^  is 
important  because,  even  if  the  pressure  response  is  ideally 
equalized  by  a  minimum  phase  network,  it  still  remains  as 
the  factor  which  causes  a  time  delay. 

The  response  is  0^taineci  from  the  Hilbert 
transform  [ 31 ] , 


00 

f,nA(u>,)du),/(u>’  -  u>)tt  . 

— 00 


(  2  .  A  5  ) 


For  the  exact  solution  of  Eq.  (2.45),  one  must  know  the 

pressure  response  A(u))  over  the  entire  region  of  u).  But  the 

range  actually  known  is  from  kA  »  0.1  to  10.0.  For  kA  * 

10.0  to  819.1,  the  pressure  responses  were  calculated  by  the 

geometrical  approximation  method.  This  frequency  region 

was  sufficient  to  obtain  the  required  precision  of  ^min^^ 

in  the  frequency  range  from  kA  =  0.1  to  10.0.  Then, 

<f>  %  (w)  is  obtained  by 

Qi.  Xp3S  S 


<t>  , .  (w) 

allpass 


♦tot.!4*0  '  'tWn<“) 


(2. 46) 
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By  dividing  Gq.  (2.46)  by  k,  one  obtains  the  corresponding 
distance  due  to  the  phase  delay,  which  will  be  referred  to 
here  as  the  "acoustic  center".  The  characteristics  of  the 
acoustic  center  for  various  values  of  H/A  are  shown  in 
Figures  2.21-2.22.  The  results  in  the  low-frequency  region 
are  not  quite  correct  because  of  using  the  discrete  Fourier 
transform  with  the  frequency  interval  of  0.1  to  calculate 
Gq .  (2.45).  Also,  a  small  error  in  the  phase  response  can 
cause  relatively  large  fluctuations  in  this  region.  The 
result  for  the  concave  dome  with  H/A  *  1.0  is  not  shown 
because  of  the  difficulty  of  taking  into  account  the  sharp 
dip  of  the  pressure  response  around  kA  -  4.15. 

It  is  interesting  to  note  that  the  position  of  the 
acoustic  center  of  a  convex  dome  is  constant  and  equal  to  H 
over  the  entire  frequency  range  from  kA  *  0.1  to  10.0.  The 
same  result  is  obtained  from  the  purely  geometrical 
approximation.  The  reason  for  the  existence  of  the  acoustic 
center  at  the  front  of  the  source  is  shown  in  Appendix  D 
for  the  case  of  a  line  source  as  an  example. 

The  position  of  the  acoustic  center  of  the  concave 
dome  is  also  constant  except  in  the  very  high-frequency 
region,  and  is  located  between  the  baffle  surface  and  the 
top  of  the  dome.  The  results  of  Figures  2.21-2.22  are  very 
important  because  they  show  the  location  of  the  source  after 
the  equalization  of  the  pressure  response.  The  equalized 
responses  of  the  convex  and  concave  domes  are  nond i s pe r s i ve 
in  the  range  used  for  an  actual  loudspeaker. 
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2.3.4  Directlvi ty  Pattern 

The  directivity  pattern  is  one  of  the  most  important 
characteristics  of  a  transducer.  If  one  knows  both  the  on- 
axis  pressure  level  and  the  directivity  pattern  at  a 
specific  frequency,  total  radiated  power  and  the  radiation 
impedance  at  that  frequency  can  be  calculated. 

Directivity  patterns  of  the  convex  dome  are  shown  in 
Figures  2.23-2.28,  along  with  those  of  a  flat  piston.  The 
on-axis  pressures  are  lower  than  the  90-degree  off-axis 
pressures  for  the  convex  domes  with  H/A  ■  1.0  or  0.75,  and 
at  kA  ■  1.0  and  2.0.  From  this  frequency  region,  the  convex 
dome  has  a  wider  directivity  pattern  than  the  flat  piston. 
The  convex  dome  has  the  narrowest  directivity  pattern 
around  kA  «  5.0  (Figure  2.26).  Again,  at  higher 
frequencies,  lower  on-axis  pressure  responses  than  the  off- 
axis  pressure  responses  are  observed  (Figures  2.27  and 
2.28). 

Figures  2.29-2.35  show  the  directivity  patterns  of  a 
concave  dome.  The  concave  dome  has  almost  the  same 
directivity  patterns  as  the  flat  piston  until  about  kA  * 
3.0.  Figure  2.32  shows  the  directivity  pattern  at  kA  = 
4.15,  which  is  the  frequency  at  which  the  on-axis  pressure 
response  has  a  sharp  dip  (Figure  2.10).  The  on-axis 
pressure  is  much  lower  than  the  66-degree  off-axis  pressure. 
This  results  suggests  that  the  normal  velocity  distribution 
on  the  opening  of  the  cavity  is  not  in  a  uniform  phase.  In 
the  higher  frequency  region,  the  directivity  pattern  of  a 
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gure  2.26  Directivity  patterns  of  a  convex  dome  an 


flat  piston  at  kA 


Figure  2.29  Directivity  patterns  of  a  concave  dome 


Figure  2.30  Directivity  patterns  of  a  concave  dome 
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Figure  2.33  Directivity  patterns  of  a  concave  dome  and 

®  _  .  it 


flat  piston  at  kA 
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Figure  2.35  Directivity  patterns  of  a  concave  dome  and 
flat  piston  at  kA  =  10.0. 


concave  dome  has  similar  properties  to  those  of  the  convex 
dome.  Neither  dome  shows  zero  pressure  at  any  angle.  As  a 
practical  point  of  view,  it  should  be  mentioned  that  the 
convex  dome  has  wider  directivity  than  the  concave  dome  or 
the  flat  piston  in  the  frequency  range  normally  used  for  the 
loudspeaker . 

2.3.5  Energy  Flow  and  P  res  sure  Distribution 

Besides  the  sound  pressure,  the  important  and 
interesting  quantity  used  to  describe  the  sound  field  is  the 
energy  flow.  Therefore,  the  sound  intensity  was  calculated 
and  the  results  are  presented  for  some  ’:A  values. 

The  x  and  z  components  of  intensity  in  the  field  are 
defined  as 

I  -  Re[p  u*  ]  /  2  .  (2.47) 

x  ,  z  X  ,  z 

Then,  the  direction  of  the  particle  motion  and  the  intensity 
IQ  in  that  direction  are  given,  respectively,  by 

6  -  tan”1  (I  /I  )  (2.48) 

x  z 

and 

?  2  1/2 

I  =  (I2  +  I2)  .  (2.49) 

o  X  z 

These  two  values  are  simply  described  by  an  arrow  with  its 
direction  and  length. 

The  graphs  of  energy  flow  and  pressure  distribution  of 
the  convex  dome  for  several  H/A  values  are  shown  in  Figures 


2.36-2.44.  The  arrow  shows  the  relative  amplitude 
proportional  to  ( 1Q ) *  ^  and  the  direction  of  IQ  at  the  point 
of  the  tail  of  the  arrow.  The  size  of  the  circle  is 
proportional  to  the  absolute  value  of  the  pressure  at  its 
center.  For  the  convex  dome  with  H/A  =  1.0  at  kA  *  1.0 
(Figure  2.42),  the  energy  tends  to  flow  from  the  high- 
pressure.  region  (near  axis)  to  the  low-pressure  region 
(off-axis)  on  the  surface  of  the  radiator.  As  a 
consequence,  at  distances  greater  than  approximately  3 
radiator  radii,  the  on-axis  sound  pressure  is  smaller  than 
the  90-degree  off-axis  sound  pressure.  As  the  frequency 
increases,  the  energy  starts  to  be  radiated  normal  to  the 
surface,  and  the  bundle  of  energy  is  confined  near  the  a::is. 

The  pressure  distribution  and  the  energy  flow  of  the 
concave  dome  are  shown  in  Figures  2.45-2.54,  The  energy 
radiated  from  the  concave  dome  with  H/A  »  1.0  at  frequency 
kA  =  1.0  (Figure  2.51)  flows  along  the  axis  inside  the 
cavity.  The  pressure  is  constant  along  the  opening  of  the 
cavity,  but  decreases  along  the  axis  even  Inside  the  cavity. 
This  shows  that  the  effect  of  the  cavity  cannot  be 
represented  by  a  simple  lumped  compliance  in  the  analog 
circuit.  At  kA  =  3.0  (Figure  2.52),  the  pressure  along  the 
opening  is  not  constant,  and  the  energy  flow  is  not  parallel 
to  the  axis  anymore.  The  graph  at  kA  =  4.15  (Figuie  2.53) 
shows  some  interesting  results.  The  energy  radiated  from 
some  area  of  the  radiator  surface  goes  around  about  360 


degrees  inside  the  cavity 


and 


goes  out  to 


semi- inf inite 


Figure  2.37  Energy  flow  and  pressure  distribution 
convex  dome  with  H/ A  =  0.5  at  kA  =  3.0. 


Figure  2.38  Energy  flow  and  pressure  distribution 
convex  dome  with  H/ A  *  0.5  at  kA  =  10.0. 


Figure  2.40  Energy  flow  and  pressure  distribution 
convex  dome  with  H/ A  =  0.75  at  kA  =  3.0. 


Figure  2.42  Energy  flow  and  pressure  distribution 
convex  dome  with  11/ A  =  1.0  at  kA  =  1.0. 


Figure  2.44  Energy  flow  and  pressure  distribution 
convex  dome  with  H/A  =  1.0  at  kA  =  10.0. 


H/A=  0.5 


Figure  2.45  Energy  flow  and  pressure  distribution 
concave  dome  with  H/A  *  0.5  at  kA  1.0. 
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Figure  2.46  Energy  flow  and  pressure  distribution 
concave  dome  with  H/A  *  0.5  at  kA  =  3.0. 


Figure  2.47  Energy  flow  and  pressure  distribution 
concave  dome  with  H/ A  *  0.5  at  kA  =  10.0 


H/A=  0.75 
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Figure  2.48  Energy  flow  and  pressure  distribution  of 
concave  dome  with  H/ A  =  0.75  at  kA  =  1.0. 


U) 

h-  CO 
®  <5> 


RADIATOR 


H/A« 


RADIATOR 


d 1st  ribut  ion 


t 


K 

■d 

k: 

A 


C  O 
o  • 

•r4  O 


9 

^  R 


« 


*9  (0 
O 


V  • 

JU 

9 

«  « 


9  <u 
o  0 
^  o 

^  *9 


0) 
>s  > 
00  <0 
U  CJ 
4>  C 
C  O 

w  o 


ir> 


0) 

h 

9 

00 


RADIATOR 


86 


space  from  near  the  rim  of  the  dome.  The  results  of  the 
soind  pressure  response  and  directivity  pattern  at  this 
frequency  are  explained  clearly  by  the  energy  flow  and 
pressure  distribution.  At  kA  •  10.0  (Figure  2.54),  the 
energy  radiated  normal  to  the  radiator  surface  converges 
around  the  center  of  the  sphere  and  then  spreads  out  to 
semi-infinite  space. 

As  was  shown  above,  the  graphical  representation  of 
the  energy  flow  as  well  as  the  pressure  distribution  is  a 
very  useful  tool  for  understanding  the  radiation  problem. 

2.3.6  Comparison  of  Results  with  other  Methods 

The  results  of  the  radiation  and  diffraction  problems 
are  compared  with  those  obtained  by  different  methods.  As 
was  mentioned  earlier,  Ikegaya  [6]  calculated  the  on-axis 
pressure  responses  of  a  convex  dome,  approximating  it  by  an 
oblate  spheroid  for  the  cases  of  H/A  ■  0.2,  0.37,  0.5,  and 
1.0.  The  last  case  corresponds  exactly  to  the  radiation 
from  the  convex  dome  with  H/A  ■  1.0.  His  results  for  H/A  * 
0.5,  and  1.0  are  plotted  in  Figure  2.55  along  with  the 
results  obtained  by  the  present  method.  Both  responses  for 
H/A  *  1.0  are  in  good  agreement  with  each  other  except  the 
small  discrepancy  in  the  low-frequency  region.  The  present 
method  Lacks  some  accuracy  in  this  region,  but  is  accurate 
enough  in  mid-  to  high-frequency  regions.  For  H/A  *  0.5, 
both  results  are  significantly  different  from  each  other  in 


the  high-frequency  region 


This  difference  is  considered 


37 


Figure  2.55  Far-field  on-axis  pressure  responses 
convex  dome  and  an  oblate  spheroid. 
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reasonable  because  the  oblate  spheroid  has  a  flatter  region 
around  the  z-axis  than  does  the  convex  dome  with  the  same 
H/A.  It  seems  that  the  oblate  spheroid  suffers  larger 
energy  flow  to  the  off-axis  direction  in  the  middle- 
frequency  region  since  it  has  smaller  normal  velocity 
distribution  around  its  rim  compared  with  the  dome. 

Figure  2.56  compares  the  on-axis  pressure  responses  of 
the  concave  dome  obtained  by  our  method  (LSM)  and  by  the 
finite  element  method  (FEM).  The  programming  and  the 
calculation  by  the  FEM  were  carried  out  by  the  staff  of  the 
Loudspeaker  Section  of  Consumer  Products  Research 
Laboratory,  Mitsubishi  Electric  Corporation  in  Japan,  at  the 
author's  request.  The  number  of  axisymmetric  triangular 
elements  were  ?.40,  360,  and  450  for  H/A  ■  0.5,  0.75,  and 
1.0,  respectively.  The  boundary  condition  on  the  opening  of 
the  cavity  was  given  in  the  form  of  self  and  mutual 
impedances  of  each  mode.  Inside  the  element,  the  pressure 
was  expressed  as  a  combination  of  pressures  at  six  nodes 
using  a  shape  function  of  quadratic  form.  The  results  of 
the  FEM  show  higher  resonance  frequencies.  But  the 
differences  of  two  curves  for  H/A  »  0.5,  0.75,  and  1.0  are 
mostly  less  than  ldB,  respectively.  The  agreement  of  these 
two  methods  indicate  that  the  present  method  is  good  enough 
for  the  discussion  of  the  radiation  characteristics  of  the 


concave  dome  radiator. 


CHAPTER  III 


DIFFRACTION  BY  CONVEX  AND  CONCAVE  DOMES 


3 . 1  Geometrical  Model 

This  chapter  deals  with  the  diffraction  of  sound  by 
the  convex  and  concave  domes  that  were  used  as  radiators  in 
Chapter  II* 

Figures  3.1  and  3.2  show  ring  sources  with  a 
concentric  convex  or  a  concentric  concave  dene  in  an 
infinite  baffle,  respectively.  The  ring  source  has  radius 
A  j  and  source  strength  2trA^Q  .  The  convex  and  concave  domes 
are  represented  as  a  portion  of  a  sphere  of  radius  R,  and 
have  height  H  and  radius  A. 

The  rectangular  coordinate  system  (x,y,z)  is  used  with 

the  origin  0.  The  polar  coordinate  system  (r,0,$)  is  used 

with  the  origin  Oj  at  the  center  of  the  sphere.  The 

surface  of  the  convex  or  the  concave  dome,  which  is  denoted 

by  S<1).  is  represented  by  (r  »  R,0<9<8  ,  0<<J><2ir)  or 

—  —  o  —  — 

(r  »  R,0  <0<r,  0<d><2iT)  ,  respectively.  In  Figure  3.2,  the 

o—  —  —  — 

opening  of  the  cavity  (  x^+y^  £  A‘,z«0)  is  denoted  by 

For  the  same  reason  as  in  Chapter  II,  the  model  shown 
in  Figure  2.1  will  be  used  for  the  diffraction  problem  by 
the  convex  dome;  however,  it  does  not  show  the  ring  v.ource. 
The  imaginary  surface  of  the  sphere  with  the  origin  Dj  is 
denoted  by  S  ^  ^ . 
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Figure  3 . 


I  Geometry  of  a  ring  source  with  a  concentric 
convex  dome  in  an  infinite  baffle. 


RING 

SOURCE 


Figure  3  . 


Geometry  of  a  ring  source  with 
concave  dome  in  an  infinite  baffle. 
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3.2  Mathematical  Discuss  ions 


3.2.1  Diffraction  by  a  Convex  Dome 

The  present  diffraction  problem  is  actually  a 
radiation  problem  from  a  ring  source  with  a  concentric 
convex  or  concave  dome  as  a  diffracting  object.  Hence,  the 
mathematical  methods  described  in  Sections  2.2.1  and  2.2.2 
are  also  applicable  to  the  present  problem. 

In  the  diffraction  problem  of  the  model  shown  in 
Figure  2.1  (with  the  ring  source  missing),  the  center  of  the 
sphere  Is  chosen  as  the  origin  of  the  coordinate  system,  and 
the  boundary  conditions  are  given  on  the  surface  of  the 
sphere  and  S^).  The  total  velocity  potential  ^(O 

outside  the  sphere  can  be  expressed  by 

«US)  "  MS)  +  <MO  >  (3- 1> 

t  s  d 

where  ijj  (£)  is  the  potential  due  to  the  ring  source,  and 
s 

^(£)  is  the  secondarily  radiated  (reflected)  potential  by 
the  convex  dome.  On  the  radiator  surface  S^),  (O  must 
satisfy  the  Neumann  boundary  condition, 


i-  £  ♦«<’••  W 


*sV[>U  ■  0  •  (3'2) 


Rewriting  Eq .  (3.2),  one  obtains 


JL  O)  1 


oil 


(3.3) 


3 
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Thus,  on  the  surface  L  ^  ^  and  f^^(0  are  given  by 


and 

f(1)(0  "  ^*s<5)  »  (3.5) 

respectively. 

The  velocity  potential  ips(r,B)  at  the  point  (r,9,0)  is 
given  by 

QA.  (2tt  -ikd 

4>s(r,e)  -  -fir  J  (e  1/d1)  d*1  ,  (3.6) 


where  d,  is  the  distance  between  the  point  of  interest 
(r,9,0),  and  the  point  of  line  element  Ajd<J>.  of  the  rin8 
source  at  (R^ ,0^ ,4>^) ,  which  is  expressed  by 


2  2  2  1/2 
d^  =  [(rsinO  -  A^cos<|»^)  +  (A^sinij^)  +  (rcos0  -  zq)  J 

(3.7) 


The  normal  velocity  distribution  on  the  surface  S 
dome  is  therefore  given  by 


(1) 


of  the 


l~  37^s(r«9>Ir«R 


QA 


1 

2 71  J 


2* 


0 


[(1  +  ikd1)e_^kdl  { (Rsin0  -  AjCOS^) 


\  sin9  +  (RcosO  -  zo)cos6}/d  j'jd^  •  (3.8) 


Expressing  the  unknown  function  i'(j(  )  in  terms  of  an 

infinite  number  of  trial  functions. 
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iK(C)  "  y  a  h  (kr)P  (cos0)  ,  (3.9) 

d  n-0  ”  °  ° 

the  function  defined  by  Eq.  (2.10)  is  obtained  from 
Eqs.  (3.4)  and  (3.9)  as 

4'(1)(6)  -  -  kh'(kR)P  (cos0),  0  <6<  0  .  (3.10) 

n  n  a  - o 

The  symmetry  of  the  sound  field  about  the  xy-plane 
requires  the  symmetry  of  the  velocity  potential  ^(£)  itself 
about  the  xy-plane;  i.e.,  the  velocity  potentials  at  point  P 
(r,0,O)  and  Pj  (r2»02»O)  are  equal  to  each  other  (potential 
matching).  From  Eq.  (3.9),  one  obtains 

oo 

7  a  (h  (kR)P  (COS0)  -  h  (kr~)P  (cos0„)}  -  0  .  (3.11) 

n-0  n  n  n  n  2  n  2 

In  the  case  of  potential  matching,  the  operator  L  is  not 
expressed  explicitly.  By  comparing  Eqs.  (2.4)  and  (3.11), 
one  comes  to  the  interpretation  that  the  operator  L  gives 
the  difference  between  two  velocity  potentials  at  a  point 
and  its  related  point  (in  this  case,  symmetric  point).  One 
can  define  ^"^(O  and  f^2^(£)  as 

tJ2)(0)  =  (hn(kR)Pn(cos0)  -  hn(kr2)Pn(cos02)]  (3.12) 

and 

f(2)(0)  =■  0  .  (3.13) 


Now,  applying  Eqs.  (3.5),  (3.8),  (3.10),  (3.12),  and  (3.13) 
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to  Eq.  (2.12),  one  obtains 


N 

l  *n[ 

n*0  n 


,<1) 


V(1)*(0)H'(l)  (0)dS  + 
m  n 


.(2) 


^2)*(0)^2)(e)dsl 

m  n 


S'*'  S 

f(1)*(0)f(1) (0)dS  ,  n  -  0,  1,  .  .  N  ,  ( 3 . 1 A  ) 


,(D 


m 


where  dS  =  2ttr  sin0dd,  and  Eq.  (2.12)  was  assumed. 


3.2.2  Diffraction  by  a_  Concave  Dome 

The  diffraction  problem  by  a  concave  dome  is  also 
treated  here  as  a  radiation  problem  from  a  ring  source  in  an 
infinite  baffle  with  a  concentric  concave  dome  (Figure  3.2). 
The  velocity  potential  inside  a  sphere  of  radius  R  with  the 
origin  at  0^  is  expressed  by 

i|>t(0  *  4JS(°  +  V0  *  (3.15) 

The  potential  i|)g(£)  is  c^e  one  radiated  by  the  ring  source 
existing  only  in  the  semi-infinite  space  on  the  right  side 
of  the  xy-plane.  The  total  potential  ^t(£)  inside  the 
sphere  can  be  represented  as 

OO 

iMO  =  l  a  j  (kr)Pn(cos0)  .  (3.16) 

n*0 

On  the  surface  of  the  concave  dome  this  must  satisfy 

the  following  condition: 


3 

3r 


'  t 


(O 


3 

3r 


V° 


0  . 


(3.17) 


This  defines  f^(Q),  and  such  as 
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(1) 


_3_ 

3r 


(3.18) 


f(1)(6)  -  0  ,  (3.19) 

and 

i(1)(6)  =  -  kj  (kR)P  (ccsQ)  .  (3.20) 

n  n  n 

In  the  region  inside  the  sphere  and  outside  the 

cavity,  the  total  potential  is  divided  into  two  components: 

the  direct  velocity  potential  ^g(£)»  and  the  diffracted 

velocity  potential  tli,(£)  due  to  the  concave  dome.  The 

d 

potential  4^  ( O  is  given  by  the  same  equation  as  Eq.  (3.6). 
8 

The  potential  1^(0  is  obtained  by  the  Rayleigh  integral  of 
the  normal  velocity  distribution  u^r,,^)  on  the  opening 
S^)f  which  is  given  by 

“n(r2'V  ■  •  1  4  Vr-8)lr  ■  r2 

e  -  e2 

-  -  I  4  »t>  If  +  M  Wc>  I  lr  -  r2  •  °-21) 

0  *  ©2 

The  derivatives  ?r/3z  and  30/3z  are  given  from  the  geometry 
by 


3r/3z  =  cosB 


(3.22) 


and 


30/3z 


-  sinO/r  . 


(3.23) 
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From  Eqs.  (3.16),  (3.22),  and  (3.23),  one  obtains 


uB(r2,62)  "  "  7  an[kcos92j^(kr2)Pn(cos92) 

n*0 


+  sin  82jn(kr2)P^(cos62)/r9]  .  (3.24) 


Then,  s^(r,8)  is  given  by 


i(j  (r,6)  -  ( 

d  j,(2) 


-ikd 

uB(r2’62)  (e  /d2)dS] /2rr  ,  (3.25) 


where  the  coordinate  of  the  surface  element  is  ^r2’^2’^2^ 
and 


2  2 

d  -  [  (rsinQ  -  r2sin92cos<|»2)  +  (r2sin02sin(j)2) 


+  (rcos0  -  r2cos02)  ] 


1/2 


(3.26) 


From  Eqs.  (3.16),  (3.24),  and  (3.25),  following  equation  can 
be  obtained: 


N 


r 2ir  -ikd. 


y  a  j  (kr)P  (cos0)  ■  QA,  [ 
nJn  n  1 

n*0 


{e  /d  1 }  deb  1 1  /  2tt 


l  a  [ 

4-  n 


n*0 


.(2) 


(kcos07j ' (kr, )P  (cos07) 
2  n  2  n 


-ikd. 


+  sin20.,j  (kr„)P' (cos07)/r7)  (e  ~/d2)dS] /2tt .  (3.27) 

2Jn  2  n  2  2  z 


The  cond i t ion 

expressed 

by  Eq.  (3.27)  must  be  satisfied 

everywhere 

in 

the 

region 

inside  the  sphere 

except  at 

the 

cavity. 

But  , 

i  f 

the  imaginary  surface  of 

the  sphere 

i  s 

chosen  as 

the 

one 

where 

the  condition  of 

Eq.  (3.27) 

i  s 
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imposed,  Che  nonuniqueness  problem  occurs  [see  Appendix  C]. 
With  this  in  mind,  it  is  chosen  that  Che  condition  of 
Eq .  (3.27)  is  required  on  the  opening  of  Che  cavity  S^2^. 
Rewrite  Eq .  (3.27)  such  as 

00 

y  a  -^^(r.e)  -  fU)(0)  ,  (3.28) 

L  n  n  n 

n-0 

where 

(r ,0)  -  1  (kr)P  (cos6)  +  [f  {kcos0  j’(kr,)P  (cos0_) 

n  n  n  J  4  n  ^  n  4 

2  -ikd2 
+  sin20J  (kr^P’Uosej/r.}  (e  z/d  )dS]/2TT  (3.29) 

4  n  4  n  44  4 


and 


f(2) 


(r,0) 


-ikd. 

(e  x/d1)dcf>1]/2TT 


(3.30) 


with  the  relationship  between  r  and  9  such  as  Eq .  (2.37). 
Now,  the  unknown  coefficients  an,  n  **  0 ,  1  ,  .  .  .,  N  are 
obtained  from  the  set  of  (N+l)  simultaneous  equations  with 
the  same  form  as  Eq .  (3.14): 


1  n  “"‘i 
n=0  4 


.(1) 


H,(1)*(0)'i'(1)(0)dS  +  f  4'(2)A(r,0)4'(2)(r,0)ds" 
in  n  J  /9 \  ^  ^ 


(2) 


4,^2^*(r,Q)f^2^(r,0)dS,  in  =  0,  1,  ,  .  ., 
m 


N  .  (3.31 ) 


3.3  Reaulta  and  Discussions 

Corresponding  to  the  radiation  problem  from  a  convex 
or  a  concave  dome,  the  he igh t /radius  ratio  of  the  dome  was 
chosen  such  that  H/A  ■  0.5,  0.75,  and  1.0.  In  most  of  the 
cases,  the  radius  of  the  ring  source  was  equal  to  1.5A. 
The  weighting  factor  q(£)  in  Eq .  (2.12)  was  kept  equal  to 
1.0  as  in  Chapter  II.  The  accuracy  of  the  ults  was 
confirmed  by  checking  the  error  factor  of  Eq .  (2.9)  at  each 
frequency.  The  maximum  value  of  index  N  was  60  for  kA  -  40. 

3.3.1  On-J xi s  Pressure  Response 

The  iar-field  on-axis  pressure  responses  of  the  ring 
source  with  the  convex  dome  are  shown  in  Figure  3. 3-3. 5, 
where  the  strength  of  the  ring  source  is  inversely 
proportional  to  u>  .  The  effect  of  the  dome  appears  from 
approximately  kA  =  1.0,  giving  gradual  rise  to  the  response 
as  the  frequency  increases.  The  height  of  the  plateau  is 
about  2.5dB  and  does  not  depend  strongly  on  the  H/A  ratio. 
The  differences  of  the  levels  of  the  peaks  and  dips  in  the 
high-frequency  region  are  more  than  lOdB  for  H/A  -  0.5  and 
0.75.  As  the  H/A  ratio  increases,  the  differences  decrease, 
and  the  peaks  and  dips  are  shifted  to  higher  frequencies. 
The  results  indicate  that  the  radiat  on  efficiency  is  highly 
frequency  dependent. 

Figures  3. 6-3. 8  and  Figures  3.9-3.11  show  the  on-axis 
pressure  responses  of  the  ring  source  with  a  concave  dome 
for  A./A  =  1.5  and  2.0,  respectively.  The  concave  dome 
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begins  Co  have  some  effect  on  the  response  from  about  the 
same  frequency  (kA  ■  1.0)  as  the  convex  dome  for  A^/A  »  1.5. 
The  first  effect,  however,  appears  as  a  large  dip  on  the 
frequency  response,  contrary  to  the  results  shown  in  Figures 
3. 3-3. 5.  The  fluctuations  of  the  response  are  much  smaller 
compared  with  those  of  the  convex  dome.  As  the  H/A  ratio 


increases , 

the 

frequencies 

o  f 

the  peaks  and 

dips  are 

lowered . 

When 

the  A  ^ 

/A 

ratio  is  increased 

(Figures 

3.9-3.11) , 

the 

effect 

of 

the 

cavity  on  the 

response 

decreases,  and  the  lowest  frequency  at  which  the  effect  is 
evident  also  decreases. 

These  results  show  that  the  object  (convex  or  concave) 
in  the  infinite  baffle  has  a  large  effect  on  the  radiation 
from  an  adjacent  source.  One  effect  of  practical  importance 
is  demonstrated  by  Figure  3.12,  which  shows  the  difference 
in  the  response  when  height  H  is  changed  by  +0.02A  from  H  * 
0.5A.  For  example,  when  A  ■  15cm,  the  change  in  the  height 
is  +  3mm.  Figure  3.12  indicates  that  the  change  of  the 
height  of  the  dome  will  produce  amplitude-modulations  of 
^bout  -30  to  -40dB  at  some  frequencies. 

3.3.2  Directivity  Pattern 

Directivity  patterns  of  the  ring  sources  are  also 
affected  by  the  convex  or  the  concave  object  in  the  infinite 
baffle  as  shown  in  Figures  3.13-3.15  and  Figures  3.16-3.18 
for  the  case  of  H/A  =  1.0,  respectively.  Figures  also 
include  some  of  the  directivity  patterns  of  the  ring  source 
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Figure  3.13  Directivity  patterns  of  a  ring  source  with 
concentric  convex  dome  and  a  flat  piston  f 
H/A  =  1.0  and  A  /A  =  1.5  at  kA  *  1.6. 


Figure  3.16  Directivity  patterns  of  a  ring  source  with  a 
concentric  concave  dome  and  a  flat  piston  for 
H/ A  -  1.0  and  A./A  *  1.5  at  kA  ■  1.2. 


I  * 


erns  of  a  ring 


117 


Figure  3.18  Directivity  patterns  of  a  ring  source  with  a 
concentric  concave  dome  and  a  flat  piston  for 
H/A  -  1.0  and  A./A  -  1.5  at  kA  -  10.0. 
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only.  Again,  the  convex  and  concave  domes  work  in  opposite 
ways.  The  convex  dome  widens  the  directivity  pattern  at  kA 
*  1.6  as  well  as  Increases  the  on-axis  response,  resulting 
in  a  larger  efficiency  (Figure  3.13).  The  directivity 
pattern  of  the  ring  source  with  the  concave  dome  at  kA  *  1.2 
is  much  narrower  than  the  one  without  the  dome  (Figure 
3.16).  In  the  high-frequency  region,  the  effects  of  the 
domes  are  too  complicated  to  be  described  simply. 

3.3.3  Ene r gy  Flow  and  Pressure  Distribution 

The  pressure  distribution  and  energy  flow  around  the 
source  at  several  frequencies  for  H/A  ■  1.0  are  shown  in 
Figures  3.19-3.24.  The  circles  indicate  the  relative 
magnitude  of  the  pressure,  and  the  arrows  show  both  the 
direction  of  the  energy  flow  and  the  square  root  of 
intensity.  The  figures  show  that  the  boundary  condition  of 
zero  velocity  distribution  normal  to  the  baffle  and  the 
surface  of  the  dome  is  satisfied.  The  results  in  Figures 
3.22-3.24  were  calculated  using  the  left-hand  side  of 
Eq.  (3.27)  inside  the  cavity,  and  the  right-hand  side  of  the 
same  equation  outside  of  it  with  the  maximum  order  of  n 
truncated  at  N . 

The  convex  dome  has  the  property  that  it  diffracts  the 
sound  more  in  the  high-frequency  region.  On  the  other  hand, 
the  concave  dome  does  not  have  much  effect  in  the  high- 
frequency  region.  At  kA  »  10,  the  energy  flows  along  the 
opening  of  the  cavity  as  if  there  were  no  cavity  (Figure 
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3.24).  An  interesting  phenomenon  is  observed  in  Figure  3.23 
which  shows  a  couple  of  vortexes  around  the  cavity  at  kA  = 
4.6.  A  similar  phenomenon  was  observed  by  Schultz  et  al.  in 
the  measurement  of  intensity  flow  in  a  reverberant  room,  and 
he  suggested  as  its  reason  a  coupling  between  different 
modes  [ 32 ] . 

Consider  a  sound  field  inside  a  room  consisting  of 
two  tangential  modes  which  is  represented  by 


*Mx,y)  =  cos(3ttx/1  )c.os(TTy/l  )  +  icos(7rx/l  )cos(2iry/l  )  . 

x  y  x  y 


(3.32) 


when  ly/lx  -  (3/8)1/2. 


1/2 

These  two  modes  have  the  same  wavenumber  (35/3)  it/1^. 

The  pressure  and  the  particle 
velocity  of  this  potential  are  obtained  by  Eqs.  (2.2)  and 
(2.3),  respectively.  Then,  the  intensity  is  given  by 
Eq.  (2.46). 


The  results  calculated  by  Eqs.  (3.32)  and  (2.46)  are 
shown  in  Figure  3.25,  which  shows  the  circulation  of  energy 
very  clearly.  In  the  case  of  a  single  mode,  there  is  no 
(time  average)  energy  flow  because  the  pressure  and  the 
particle  velocity  are  90  degrees  out-of-phase  with  each 
other.  If  the  field  is  of  multiple  modes,  however,  energy 
flows  from  one  mode  to  another,  causing  the  circulation  of 
energy.  If  a  semi-infinite  space  contains  any  part  which  is 
resonant,  a  similar  phenomenon  occurs.  In  the  present  case, 
the  cavity  behaves  as  a  resonant  subspace,  producing  a 
couple  of  vortexes  of  energy  flow  at  specific  frequencies. 


CHAPTER  IV 


CONCLUSIONS 


4.1  Remarks  on  the  Mathematical  Method 
The  least  square  error  method  is  found  to  be  quite 
successful  for  solving  the  radiation  and  diffraction 
problems.  The  flexibility  of  this  method  allows  it  to 
handle  the  complex  shape  of  a  diaphragm  or  a  diffracting 
object.  This  was  shown  in  the  radiation  and  the  diffraction 
problems  by  a  concave  dome,  where  the  boundary  condition  was 
given  on  the  dome  surface  and  on  the  opening  of  the  cavity. 
Following  are  the  interesting  techniques  employed  in  the 
application  of  this  method  to  the  preceding  problems,  which 
worked  quite  well  in  improving  the  accuracy  of  the  results: 

1.  The  most  appropriate  point  was  chosen  as  the  origin  of 
the  coordinates. 

2.  The  boundary  conditions  were  given  on  the  imaginary 
surface  as  well  as  on  the  radiator  surface. 

3.  New  types  of  boundary  conditions  were  used,  which  were 
systematically  treated  in  the  same  way  as  the  Dirich let 
and  Neumann  boundary  conditions. 
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The  counterpart  of  this  method  Is  the  finite  element 
method.  The  advantages  and  disadvantages  of  the  two 
methods  are: 

1.  In  this  method,  the  sound  field  is  expressed  as  a  finite 

number  of  terms,  each  of  which  satisfies  the  wave 

equation.  Then,  the  error  is  concentrated  on  the 
boundary  surface.  This  may  simplfy  the  evaluation  of 
errors.  In  the  finite  element  method,  sources  of  error 
exist  both  on  the  boundary  surface  and  in  the  field. 

2.  In  this  method,  one  must  be  aware  of  the  nonuniqueness 

problem  when  he  deals  with  a  sound  field  including  the 

origin.  This  must  be  avoided  by  a  proper  choice  of  the 
boundary  surface.  Normally,  the  finite  element  method 
does  not  have  this  kind  of  problem  for  radiation  and 
aif fraction  problems. 

3.  This  method  is  closer  to  the  analytical  method,  giving 
better  understanding  of  the  phenomenon. 

4.  The  largest  deficiency  of  this  method  lies  in  the  nature 

of  the  spherical  Bessel  and  Neumann  functions.  One  must 
-  1  v  very  large  and  very  9mall  numbers  at  the  same 

time  for  small  argument  of  these  functions.  The  finite 
element  method  usually  does  not  have  this  kind  of 


problem • 
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4.2  Summary  of  Results 

The  main  purpose  of  this  study  was  to  Investigate  the 
difference  of  radiation  and  diffraction  phenomena  due  to  the 
difference  of  the  diaphragm  shape.  Using  the  least  square 
error  method,  several  important  results  were  obtained,  which 
may  be  summarized  as  follows: 

1.  The  convex  and  concave  domes  have  contrasting  radiation 
and  diffraction  phenomena. 

2.  The  convex  dome  becomes  more  inefficient  as  a  radiator 
as  H/A  increases.  The  concave  dome  is  more  efficient  for 
higher  H/A  in  the  range  normally  used  for  a  loudspeaker. 

3.  The  geometrical  approximation  method  gives  a  good 
estimation  of  the  response  of  the  convex  dome  in  the 
high-frequency  region.  It  is  a  rough  estimation  for  the 
concave  dome  in  the  same  region. 

4.  The  radiation  resistance  of  the  convex  dome  is  quite 
smooth.  The  radiation  resistance  of  the  concave  dome 
has  a  large  peak  at  its  first  cavity  resonance. 

5.  The  radiation  mass  of  the  concave  dome  is  equal  to  the 
sum  of  the  mass  of  the  fluid  inside  the  cavity  and  the 
radiation  mass  of  the  flat  piston  with  the  same  radius. 

6.  The  radiation  resistance  curve  closest  to  the  +6dB/oct. 
line  is  obtained  by  a  concave  dome  with  H/A  less  than 
0.5. 

7.  The  location  of  the  sound  after  the  equalization  of  the 
pressure  response  exists  at  the  top  of  the  convex  dome 
and  slightly  behind  the  opening  of  the  cavity  of  the 
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concave  dome. 

8.  The  convex  dome  has  a  wider  directivity  pattern  than  the 

flat  piston  or  the  concave  dome.  Neither  of  them  has 

any  zero  pressure  angle,  and  they  have  similar 
directivity  patterns  in  the  high-frequency  region. 

9.  The  representation  of  sound  field  by  the  energy  flow  was 
found  to  be  very  powerful  in  understanding  the  radiation 
and  diffraction  phenomena. 

10.  The  convex  dome  diffracts  the  sound  of  the  concentric 
ring  source  much  more  than  the  concave  done  does. 

11.  The  results  indicate  that  a  cone-type  woofer  may 
diffract  the  sound  emitted  from  a  tweeter  sufficiently 
to  produce  a  discernible  amplitude-modulated  distortion. 

12.  In  the  low-  to  mid-frequency  region,  the  convex  dome 
widens  the  directivity  pattern  of  the  ring  source, 
whereas  the  concave  dome  works  in  the  opposite  way. 

13.  The  existence  of  vortexes  of  the  energy  flow  In  the 
resonant  subspace  was  observed,  and  it  was  explained  as 
the  result  of  interaction  between  different  modes. 
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4 . 3  Future  Work 

The  results  obtained  so  far  are  rather  basic  knowledge 
about  the  radiation  and  diffraction  phenomena.  For  further 
information  for  the  improvements  of  a  loudspeaker,  more 
practical  aspects  of  the  loudspeaker  should  be  taken  into 
account  in  the  calculation.  The  following  is  recommended  as 
future  work. 

4.3.1  Moda 1  Vibration  of  the  Diaphragm 

The  high  end  of  the  reproducing  range  of  an  actual 
loudspeaker  is  limited,  mainly  by  the  modal  vibration  of  the 
diaphragm  rather  than  the  acoustical  reason.  The  assumption 
of  piston-like  motion  of  the  diaphragm  is  then  not  so 
preferable  to  discuss  the  responses  in  the  high-frequency 
region  of  the  loudspeaker.  The  velocity  distribution  of  the 
diaphragm  at  its  modal  vibration  can  be  obtained  either  by  a 
measurement  or  by  a  theoretical  calculation.  Applying  this 
velocity  distribution  into  Eq.  (2.6),  the  modal  vibration  is 
easily  taken  into  account  in  the  present  method. 

4.3.2  Diaphragm  Shape  other  than  j*  Dome 

A  cone-type  loudspeaker  usually  has  a  straight  cone 
with  a  dome-type  dust  cap  at  its  center.  In  this  case,  the 
concave  dome  is  not  really  a  good  approximation  of  the  cone- 
type  loudspeaker  diaphragm.  In  order  to  investigate  the 
radiation  or  the  diffraction  by  the  cone  diaphragm,  an  exact 
shape  must  be  considered,  especially  in  the  high-frequency 
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region.  As  was  mentioned  in  Section  4.1,  the  present  method 
does  not  necessarily  depend  on  the  dome  shape.  When  the 
present  method  is  applied  to  the  diaphragm  shape  other  than 
the  convex  or  the  concave  dome,  it  is  recommended  to  choose 
the  center  of  a  sphere  as  the  origin  of  the  coordinate, 
where  the  portion  of  the  sphere  represents  the  shape  of  the 
diaphragm  most  effectively. 

Calculations  of  the  characteristics  of  some  kinds  of 
diaphragm  shapes  will  give  the  general  knowledge  about  the 
relation  between  the  diaphragm  shape  and  several  important 
responses  such  as  the  on-axis  pressure  response,  radiation 
impedance,  and  directivity  pattern.  The  interpretation  of 
this  relation  will  yield  information  about  the  "best"  shape 
of  the  loudspeaker  diaphragm,  which  may  be  different  for 
different  purposes.  It  is  needless  to  say  that,  in  the 
process  of  the  evaluation  of  the  characteristics,  the 
psychological  aspects  of  them  must  be  taken  into  account. 

4.3.3  Three-Dimensional  Problem 

The  vibration  of  a  diaphragm  with  radial  modes  is 
considered  less  effective  on  most  of  the  radiation 
characteristics.  Hence,  the  assumption  of  the  axisymmetry 
of  the  sound  field  seems  reasonable.  When  one  discusses  the 
diffraction  problem,  however,  he  should  be  able  to  deal  with 
the  problem  as  a  three-dimensional  case.  The  process  to 
extend  the  present  two-dimensional  program  to  a  three- 
dimensional  program  is  rather  a  straightforward  task.  The 
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problems  one  may  encounter  are  related  to  the  limitation  of 
the  memory  site  of  the  central  processing  unit  of  a  computer 
and  the  computation  time.  These  limitations  may  be  overcome 
by  a  large  computer  system  and  improved  algorithms  that  are 
now  available.  The  results  will  bring  some  more  important 
knowledge  about  the  diffraction  phenomenon. 


A 
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m  ■  0,  1 ,  .  .  . ,  N  . 


(A. 5) 
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APPENDIX  B 


ANOTHER  SOLUTION  TO  RADIATION 
PROBLEM  FROM  A  CONVEX  DOME 


Another  solution  rather  than  the  one  discussed  in 
Section  2.2.3  is  possible.  The  details  of  this  method  and 
its  deficiencies  are  discussed  here. 

The  sound  field  is  expressed  in  terms  of  the 
spherical  Hankel  functions  and  Legendre  functions  with  the 
origin  at  0  in  Figure  2.1: 

N 

iKr'.e')  =  l  a2nh2n(kr')P2n(coS0,)  ’  (B.l) 

n=0 

where  only  the  even  orders  of  terms  are  used  because  of  the 

symmetry  of  the  sound  field.  The  normal  velocity 

distribution  u_(8)  is  obtained  by  differentiating  Eq.  (B.l) 
n 

in  terms  of  r  such  as 


u  (0)  =  -  •  Ill  +  -ii_  .  lii) 

— V  ^  VaT-'  T  5Q'  '  ' 


u  9r  3r  30'  3r 

The  relationship  among  r',  0?,  r,  and  0  are  given  by 


(  B .  2  ) 


,2  2,2  _ 

.  =  r  +  z  2rz  cos0 

o  o 


(  B  .  3  ) 


and 


(rcosO  -  zo)tar.9'  =  rsin8  , 


(B.4) 


where 


From  Eqs.  (B.3)  and  (B.4),  one  obtains 


*  (r  -  zocos0)/r'  -  cos(0*  -  0) 

and 

30 1 

“  -  sin(0'  -  0)/r*  . 

Substitution  of  Eqs.  (B.6)  and  (B.7)  into  Eq.  (B.2) 

N 

u_(0)  =  -  [  a  [kh,  (kr')P.  (cos0,)cos(e'  -  0) 

n  n=0  n  zn 

+  h2n(kr,)P2n(cose,)sin0,sln(0;  '  e>/r * ]  . 

N°w,  f2n(e)  and  f(0)  are  defined  such  as 

*2n<0)  =  *  tkh2n(kr-)P2n(cos0‘) 

+  h_  (kr')P'  (cosO'JsinO'sinC©'  -  0)/r'] 
zn  Zn 

and 


f (0)  =  U  cos0. 

o 


Then,  the  unknown  coefficients  an,  n  -  0,  1,  . 

obtained  by  solving  (N+l)  simultaneous  equations: 


N 


l 

n=0 


a2nl 


°Y*  (0)Y_  (0)27TR2sin(0)d0] 
zm  zn 

0 

■®0  A  T 

^2m^0^UoCOS  (®)2ttR  sin(0)d0 

0 

m  =  0,  1,  2,  ....  N 


(B.6) 

(B.7) 
gi  ves 

(B.8) 

(  B .  9  ) 

(B.  10) 

,  N  are 


(B.ll) 
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This  method  looks  simple  because  it  does  not  require 
the  introduction  of  the  imaginary  surface*  But  it  was  found 
that  this  method  has  two  basic  problems.  Firstly,  the 
convergence  of  Eq .  (B.8)  to  Eq.  (B.10)  is  very  slow.  This 
is  because  the  origin  0  is  not  suitable  for  the  expansion  of 
sound  field  in  the  form  of  Eq .  (B.l).  For  the  radiation 
problem  from  a  portion  of  a  sphere,  it  seems  reasonable  to 
choose  a  coordinate  system  with  its  origin  at  the  center  of 
the  sphere.  Secondly,  the  argument  kr '  becomes  small  near 
the  z-axis  for  small  H/A  ratio.  The  real  and  imaginary 
parts  of  the  spherical  Hankel  function  become  very  small 
and  very  large,  respectively,  for  small  argument  when  the 
order  n  become  large.  This  causes  the  underflow  and 
o ve r f low  problem  in  the  process  of  numerical  computation. 
This  nature  of  the  spherical  Hankel  function  prohibits  from 
taking  higher  order  of  n  when  one  tries  to  have  good 
approximation  of  Eq .  (B.8)  to  Eq  .  (B.ll). 


‘  -!r 
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APPENDIX  C 


NONEXISTENCE  AND  NONUNIQUENESS  PROBLEMS 
OF  A  SERIES  EXPANSION  OF  THE  SOUND  FIELD 


i.  Radiation  F rom  £  Sphere 

The  nonexistence  and  nonuniqueness  properties  of  the 
solution  of  the  sound  radiation  problems  were  discussed 
theoretically  by  Copley  [11]  and  Schenck  [12].  These 
properties  are  discussed  here  focussing  especially  on  the 
series  expansion  of  the  sound  field  inside  and  outside  the 
sphere . 

The  velocity  potential  outside  the  sphere  with  only  an 
outgoing  axisymmetric  wave  can  be  expressed  by 


N 

'Kr.Q)  =  y  ah  (kr)P  (cos0)  .  (C.l) 

T  ‘  0  n  n  n 

n=0 

Generally,  the  boundary  condition  is  given  either  by  the 
pressure  distribution  or  by  the  normal  velocity  distribution 
such  as  P(0)  or  U(9)  on  the  surface  of  the  sphere. 

Then  , 

N 

p(e)  =  y  p  .  p  (cose)  ( c . 2 ) 

n  n  n 

n=0 

a  nd 


U(6) 


where 


N 

y  u  .  p  (cose)  , 

n  n 


P 


n 


(n  +  1/2) 


■7T 

P(6)P  (cosB)  sin(0)d6 
n 

0 


(  C  .  3  ) 


(  C  .  4  ) 
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a  -  -u  /kh'(kR)  .  i  ( C  .  7 ) 

n  n  n  , 

i 

Now,  h  (x)  and  h'(x)  can  never  be  zero  for  x  ^  0, 
**  n 

because  the  real  part  jn  (  x  )  and  imaginary  part  nn  (  x)  °f 
hn(x),  and  their  derivatives  satisfy  the  following  equation 
[33]  : 

j  (x)  •  n’(x)  -  n  (x)  •  j  *  (x)  =  1/x"  .  (C.8) 

n  n  n  n 

Thus,  an  is  uniquely  determined  for  any  prescribed  pressure 
distribution  or  normal  velocity  distribution. 

But,  when  the  velocity  potential  is  expressed  such  as 
(corresponding  to  the  solution  by  Green's  function  [34]) 

N 

iKr,9)  =  J.  anj  (kR)h  (kr)P  (cos9)  ,  (C.9) 

n=0  n  n  n 

the  unknown  coefficient  an  is  given  by 

a  =  p  /jwOj  (kR)h  (kR) 
n  n  n  n 


(C.  10) 
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APPENDIX  D 


ACOUSTIC  CENTER  OF  A  LINE  SOURCE 


The  acoustic  center  of  a  line  source  which  is  lying  in 
the  z  direction  between  z  =  0  and  z  -  iQ  is  calculated  as 
f  ol lows . 

The  far-field  on-axis  pressure  response  of  the  source, 
which  is  normalized  to  a  point  source  with  the  same  source 
strength  at  z  »  0,  is  given  by 


1 


-  [sinkl  +  i(l  -  coski  )]/kl 
o  o  o 

■  sin(ft)elSVft  ,  (D.l) 


where  ft  =  klQ/2. 

The  imaginary  part  of  a  minimum  phase  network  having 
the  real  part  of  Eq.  (D.l)  is  obtained  using  the  Hilbert 
transform  such  as  (using  the  integral  formula  in  [35]) 


Im(pmin) 


sin2ft'  dft' 

2ft'  ir(ft'  -  ft) 


=  (cos2ft  -  l)/2ft  ,  (D .  2 ) 

Then,  the  total  response  of  the  minimum  phase  network  is 
given  by 

Pmin  =  sin(ft)e“in/ft. 


(D.3) 


JA 


2 


h 


■*  i 
•  *1 


g 

H 


Li 


The  phase  responsv  <t>allpass  defined  by  Eq.  (2.46)  is  given 
by  the  difference  of  phase  responses  of  Eq .  (D.l)  and 
Eq .  (D.3)  such  as 


^allpass^ 


-  0  -  (-ft)  -  kl 


(D.4) 


Dividing  Eq.  (D.4)  by  k.  one  obtains  the  acoustic  center  1Q 
which  is  constant  with  ro.spect  to  frequency. 
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